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We search for electron anti-neutrinos (z7e) from long and short-duration 
gamma-ray bursts (GRBs) using data taken by the KamLAND detector from 
August 2002 to June 2013. No statistically signihcant excess over the background 
level is found. We place the tightest upper limits on W fluence from GRBs be¬ 
low 7 MeV and place hrst constraints on the relation between W luminosity and 
effective temperature. 

Subject headings: Neutrino, GRB, KamLAND 


1. Introduction 

Gamma-ray bursts (GRBs) are the most luminous phenomena in the universe. The 
duration of GRBs (AtcRs) varies in the range between 10 ms and 1000 s, with a roughly bi- 
modal distribution for so-called long GRBs of AfcRS ^ 2 s and short GRBs of AtcRB ^ 2 s. 
The progenitors of most short GRBs are widel y thought to be mergers of neutron star- 
neutron star or black hole-neutron star binaries (iMeszarosI 120061) . A favored mod el of long 
GRB progenitors is a catastrophic collapse of a massive star into a black hole flMeszaros 


20061 ). These models are supported by observa t ions of af 


galaxies for short GRBs fjVillasenor et ah 


novae associated with long GRBs fjWoosley et ah 


2005 


Fox et ah 


1999 


erglow s and identihcation of host 


20051). and obser vations of super- 


Hjorth et al.l 120031) . Both scenarios 


would result in the formation of a compact rotating black hole with an accretion disk at 
MeV or higher temperatures, which generates collimated relativistic hreball jets leading to 
GRBs. Although such a hreball model is promising and attractive, the initial condition and 
generation mechanism of the hreball jets are still unknown, since it is difficult to observe the 
optically thick center region of GRBs by electromagnetic waves. 

A potential scheme to directly expl ore the GRB center reg i on is the use of therm al neu¬ 


trinos and gravitational waves (GWs) (jSuwa fc Murasell2009t ISekiguchi et ahiDOllI) . since 


they have strong transmissivity. Thermal neutrinos are sensitive to thermodynamic prohles 
of the accretion disk, and GWs are sensitive to the dynamics of progenitors. Both are comple¬ 
mentary observations to probe GRBs. Super-Kamiokande (SK) and Sudbury Neutrino Ob- 


upper huence limits 

(Fukuda et al. 

2002: 

Aharmim et al 

201^ 

1). 0 

;hers placed limits on high 

enerev neutrinos produced in the hreball iets (Achterbere et ah 

2008: Thrane et a . 

2009: 

Abbasi et al. 2011: 

Vieregg et al. 

2011; 

Avrorin et al. 

2011 

Adrian-’ 

vlartinez et al. 

2013). 

GWs from GRBs were studied by a GW detector network {I 

\.asi et al. 

2014). 


In this paper, we present a study of electron anti-neutrinos (z/g) of a few tens of MeV 
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energy pr oduced by thermal processes from the GRB center region, especially the accre¬ 
tion disk ( Naeataki fc Kohri 2002 : Sekiguchi fc Shibata 2011 : Caballero et al. 2009) with 
the Kamioka Liquid Scintillator Anti-Neutrino Detector (KamLAND). We constrain the 
relation between the luminosity (L) and effective temperature (T) as well as W ffu- 
ence for the first time. The L-T relationship can be used to directly compare with the¬ 
oretical predictions. These limits and constraints are established using redshift-measured 
GRBs. We adopt the s tandard ACDM cosmology with Dm = 0.315, Da = 0.685, and 
Hq = 67.3 kms“^ Mpc“^ flAde et al.l 120141) throughout this paper. 


2. KamLAND detector 

The KamLAND detector is located ~lkm under the peak of Mt. Ikenoyama (36.42°N, 
137.31°E) near Kamioka, Japan. The 2,700 meter-water-equivalent (mwe) of vertical rock 
overburden reduces the cosmic-ray muon flux by almost five orders of magnitude. A schematic 
diagram of KamLAND is shown in Figured] The primary target volume consists of Ikton 
of ultra-pure liquid scintillator (LS) contained in a 13-m-diameter spherical balloon made of 
135-/im-thick transparent nylon ethylene vinyl alcohol copolymer (EVOH) composite film. 
The LS consists of 80% dodecane and 20% pseudocumene (1,2,4-trimethylbenzene) by vol¬ 
ume, and 1.36 ± 0.03 g 1~^ of the ffuor PPO (2,5-diphenyloxazole). A buffer comprising 
57% isoparaffin and 43% dodecane oils by volume, which fills the region between the balloon 
and the surrounding 18-m-diameter spherical stainless-steel outer vessel, shields the LS from 
external radiation. The specific gravity of the buffer oil is adjusted to be 0.04% lower than 
that of the LS. An array of photomultiplier tubes (PMTs)—1,325 specially developed fast 
PMTs masked to 17-inch diameter and 554 older 20-inch diameter PMTs reused from the 
Kamiokande experiment flKnme et aLlll983l) — are mounted on the inner surface of the outer 
vessel, providing 34% photocathode coverage. This inner detector is shielded by a 3.2kton 
water-Cherenkov veto detector. 

KamLAND uses the inverse beta-decay reaction to detect W: 


z/g -f- p ^ -|- n. 


( 1 ) 


This process has a delayed-coincidence (DC) event-pair signature which offers powerful 
background suppression. The energy deposited by the positron, which generates the DC 
pair’s prompt event, is the sum of the e"*" kinetic energy and annihilation 7 energies, Ep{= 
Te+ + 2me), and related to the incident W energy by = {Ee+) -|- J -|- E^^E^^/rUp, 
where E^^ and E^^ ar e neutrino and electron energy in the center of mass frame, and 


6 = {i7i^—rrip—rril)/2rnp fjStrumia &: Vissanill2003l) . In the low energy < 20MeV) range. 
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we can approximate the above relation by = Ep + 6E, where 6E = 0.78 MeV. We use this 
approximation also above 20 MeV and comment on the associated uncertainty later. The 
delayed event in the DC pair is generated by a 2.2 (4.9) MeV y-ray produced whe n the neu¬ 
tron captures on a proton (^^C). The mean neutron capture time is 207.5 ±2.8 /is (lAbe et ah 


20101 ). The angular distribution of the positron emission is nearly isotropic. Unlike in a water 
Cherenkov detector, the scintillation light is also isotropic. As a result, the positron signal 
does not provide the incoming W source direction. Due to the extremely low cross section of 
Vg, the Earth does not shadow MeV-energy extraterrestrial W- The detector therefore has 
isotropic sensitivity to GRBs. 

The event energy and vertex reconstruction are based on the timing and charge distri¬ 
butions of scintillation photons recorded by the PMTs. The reconstruction algorithms are 
calibrated with on-axis and off-axis radioactive sources deployed from a glove box installed 
at the top of the detector. The radioactive sources are ^°C o, ^^Ge, ^^Zn, ^"^^Am^Be, 


137 c 


Cs, and ^^^Po^'^C, providing energy and vertex calibration (iBerger et ah 


2009: Banks et ah 


20151 ). The overall vertex reconstruction resolution is ~ 12 cm/y^A(MeV) and energy res¬ 
olution is 6.4%/'^i7(MeV). The energy reconstruction of positrons with Ep > 7.5 MeV 
> 8.3 MeV) is verihed by using tagged ^^B (3 -decays generated via muon spalla- 


i.e. 


tion ( Abe et al.ll2010l) 


In Se ptember 2011 , the K amLAND-Zen double-beta (/5/5) decay search experiment was 
launched flGando et al.N2012l) . This experiment makes use of KamLAND’s extremely low 
background. The KamLAND detector was modihed to include a (3{3 source, 13 tons of Xe- 
loaded liquid scintillator (Xe-LS) contained in a 3.08-m-diameter inner balloon (IB), at the 
center of the detector. 


Event selection 


3.1. KamLAND DC events 


In this analysis, we use KamLAND data collected from August 3, 2002 to June 4, 2013. 
During the majority of this period, K amLAND was measuring W from nuclear power plants 


with a spectrum up t o about 8 MeV flGando et al.ll2011al. 


Earth’s deep interior (lAraki et al.ll2005l: iGando et al.ll2011b 


20131) and geological Ug from the 


20131) . Following the Fukushima 


reactor accident in March 2011, all Japanese reactor were subject to a protracted shutdown. 
The data-set is divided into two periods. Period I refers to data that was taken until the IB 
installation in September 2011. Period n refers to the data taken after the IB installation, 
which mostly coincided with the low reactor Vg flux. 
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In Period I, we search only for Vf. events with 7.5 MeV) <E^< 100 MeV which 

corresponds to the energy range of interest for GRBs with almost zero contamination from 
the reactor Vf. flnx. Dnring Period H, the reactor signal is minimal, allowing a redaction of 
the energy threshold to = 0.9 MeV. 


For the DC event pair selection, we apply the following series of selection cnts: the 
prompt energy is required to be < Ep < 100 MeV in Period fc, and the delayed energy 
to be 1.8 MeV < E^ < 2.6 MeV for neutron capture on protons or 4.4 MeV < i7d < 5.6 MeV 
for neutron capture on a hducial volume cut of i? < 6 m from the center of the balloon 
on both prompt and delayed events, a spatial correlation cut of Ai7 < 1.6 m and a time sep¬ 
aration cut of 0.5 /is < AT < 1.0 ms. Spallation cuts were used to reduce backgrounds from 
long-lived isotopes, e.g., ®Li (r = 257 ms and Q = 13.6 MeV), that are generated by cosmic 
muons passing through the scintillator. In Period II, we have to use an additional spatial 
cut for delay ed events to avoid backgrounds from the IB and its support material as shown 
in Figured] (jCando et al.l 120131) and a second-level cut using a likelihood d iscriminator to 
reduce accidental backgrounds in the low-energy region flGando et al.ll2011all . The selection 
efficiency (e^) is evaluated from Monte Garlo simulation separately for Period I {k=l) and 
Period II (/c=n) due to these additional cuts. Note that e® depends on Ep because of the 
energy-dependent second-level cut. The number of target protons in 77 < 6 m is estimated 
to be At = (5.98 ± 0.12) x 10=^F 


The total livetime during Period I was 6.91 yr and 55 DG events were observed during 
this period. In Period II, KamLAND found 88 DG events with 1.2 yr livetime. The livetime 
is dehned as the integrated period of time that the detector was sensitive to Vf. and includes 
corrections for calibration periods, detector maintenance, daily run switch, etc. The event 
rates are 9.1 x 10““^ and 8.4 x 10“^ events per hour in Period I and n, respectively. 


3.2. GRB events 

We use GRB events observed by one or more of SWIFT, HETE-2, Ulysses, INTEGRAL, 
AGILE, MAXI, and T7777M7 based on The Gamma-ray Goordinates Networl|^. Initial se¬ 
lection criteria are the requirement that the GRB be in the time period between August 3, 
2002 and June 4, 2013 and the existence of redshift and GRB-duration time measurements. 
At this stage, 256 long GRB and 21 short GRB events are left. Subsequently, all the Kam- 


^http: //gen. gsf c . nasa. gov/ 
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LAND run^ that include GRB events must have passed basic quality criteria (e.g., not a 
calibration run and stable operation). This leaves 175 long GRBs and 17 short GRBs in 
Period 1. Period n contains 38 long GRBs and one short GRB. One can see our GRB list by 
the websit^. 


4. Data analysis 

The average number of DG events and GRB events per three months is shown in Figure 
[2J In this hgure, one can see a “step” before and after the launch of the SWIFT satel¬ 
lite (November 2004) for GRB events. In contrast, there is no time dependence of the DG 
event rate during each periods. We therefore decided to analyze the whole KamLAND data 
regardless of the GRB event rate. 


4.1. Coincidence analysis 


We conduct a time-coincidence analysis between the redshift-measured GRB samples 
and the KamLAND DG events for long and short GRBs. The coincidence search time- 
window between a GRB event and a KamLAND DG event is dehned as: —fp -|- Tgrb < 
Tdc < Tqrb + ^^GRB + fp + where Tdc and Tqrb are the absolute times of the 

KamLAND DG and GRB events, respectively. A^grb is the measured GRB duration time, 
fp is 150 sec corresponding to a model-dependent, but reasonable time difference between 


the therma l neutrino production and the GRB photon production fjSekiguchi fc Suwal 12012 


Tomall2014j) . ar id F(z) is the relativistic flight-time d elay of MeV neutrinos due to non-zero 


neutrino mass flLi et aLll2005l: iGhoubev fc Kindl2003l) : 


U{z) = 


1 ml 


dz' 


2 Ei 


Iq (1 -f Z^)^i7o'\/^A + (1 + R)^Dr 


( 2 ) 


with the assumption of = ^heaviest = 87.2 meV from '^m^, < 0.23 eV (lAde et al.l 1201411 


and > 8.3 MeV in Period I and E-p^ > 1.8 MeV in Period II. All parameters of the 
time-window are hxed before the coincidence search. The total window length for the long 
GRBs is 25.2 hours (18.3 hours in Period I and 6.82 hours in Period n). The short GRBs 


^ KamLAND data-taking is stopped and restarted every day to ensure smooth data taking. The length 
of a KamLAND run is typically 24 hours long. 

^see http://www.awa.tohoku.ac.jp/KamLAND/GRB/2015 
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sum to a total of 1.45 hours of on-time window (1.33 hours in Period I and 0.11 hours in 
Period H). 

No coincidence DC events were found in the above time window for both long and short 
GRBs. We estimate the expected accidental coincidence of DC events to be 7.4 x 10“^ and 
2.2 X 10“^ for long and short GRBs, respectively. For long GBRs, the background spectrum 
is shown in Figure[3]with several expected spectra from our 90% upper limits (see l4.3p . In the 
absence of a signal, the Feldman-Cousins upper (FC) limits on the DC events are Nqq = 2.365 
and 2.435 with 90% conhdence level (CL) for long and short GRBs, respectively. 

If we use a much longer, exotic, time window, e.g., tp = 6h, four coincidence DC events 
are found for long GRBs. However, the expected accidental coincidence of DC events is 3.4. 
There is therefore no statistical evidence for the detection of W from long GRBs. 


4.2. Fluence upper limits 


There is no established neutrino production model for GRBs. We translate our FC 
limits to model-independent upper limits on W fluence, at the detector using a 

Green’s function, which represents the upper limits on monoenergetic neu trinos at that 


speci hc energy. We use the same me thodology to estimate as SK flFukuda et al. 

2002 ) and SNO f Aharmim et al. 2014 ): 


'i{Ev.) = 


N< 


90 




( 3 ) 


where is the number of GRBs and Ik{EpJ is the effective nnmber of DC events per 

one GRB with a monoenergetic spectrum in the period k\ 


"lOOMeV 


h{Ej,J = iVT 
and 


'Et 

low 




^(™xp,™s)^ 


\/2^ct(E“P) 


exp 


(^exp _ ^vis)2 


2a2(F;“P) 


( 5 ) 


et is the mean livetime-to-runtime ratiqj and E^^'^ and 77”® are the expected and measured 
prompt energies, respectively. criBD(-E) is the differential cross section of the inverse beta 
decay. cr{E) corresponds to the energy resolntion of 6.4%/A/i7(MeV). 


^Runtime is the total time of data taking. 














The 90% CL upper limits on from Kam LAND are shown fo r both lon g and short 


GRB s together in Figure |4] with the results from SK (iFnknda et alJl2nn2h and SNO (lAharmim et ah 


2014J) . Note, the results from SK and SNO treated long and short GRBs as the same. Below 


7 MeV, our analysis provides the best limits so far. 


4.3. Constraint on Inminosity and effective temperature (L-T) 


NgQ can be translated to constrain the V^s luminosity (L) and effective temperature 
(T) in the accretion disk using the assumption that the Fg flux follows the Fermi-Dirac 
distribution described: 




120 L 


El 


exp{EpjT) + 1 
The expected total flux at the detector is in Period k, 


4'‘(B„„r,L)= 


i£k 


1 + Zi 

AttS 


' 0((1 + Zi)Ey^,T,L), 


( 6 ) 


( 7 ) 


where Zi and di are the redshift and luminosity distance of the fth GRB. The luminosity 
and effective temperature upper limits (Tup, Tup) are then connected to Nqq\ 


k 


nOO MeV 


'Et 

low 


/'(Tup,Tup,Tf)dF™, 


( 8 ) 


where /(, is the visible spectrum of the DG events: 


/UUp.Lup.Bp™) = 


-100 MeV 


'Et 

low 


iVT€,£‘(i5p™)ffiBD(£i>.)'l'‘(Bp., r„p, L„p)fl(B“-, 


( 9 ) 

With the assumption of E-p^ = Ep + SE, the results obtained from KamLAND are shown in 
Figure [5l The upper limit spectra (Yl Tup, T™)) with Tup = 5,10,15 MeV are shown 

in Figure [3l 

The limits are six orders of magnitude higher than the supe rnovae W luminosi t y and 
several orders of magnitude higher than theoretical studies predict. iNaeataki fc Kohril (120021) 
analytically show that a collapsar emits W with T = 10^^ erg and T = 5 MeV in a total 
accretion mass o f 30M f7i, a initial mass of SMq, and a mass accretion rate of O.IMq/s. 


Gaballero et al.l (120091) numerically predict T = 3.5 x 10^^ erg during 0.15 sec with T = 


7.5 MeV for black hole-neutron star mergers. Here, we assumed the averaged Uf, energy 
corresponds to 3.15T. Recently, Sekiguchi presented T = 1.5-3 x 10^^ erg/s during 2-3 sec 
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with an averaged v^. energy of 20 -30 MeV for a merger of b inary neutron stars using state- 
of-the-art numerical simulations (jSekignchi fc Shibatall201lh . 


Finally, we comment about the approximation, E-g^ = + 6E. Above 20 MeV, this 

approximation is not suitable. In addition, the effect of the recoiling neutron (£'„) to Ep is 
no longer negligible. This effect adds a substantial energy bias, ~ 10%, but the uncertainty 
of Lup is much smaller than 10%. The amount of the error has no impact on our result and 
discussion. 


5. Summary 


We hnd no evidence for W associated with our sample of GRBs in KamLAND. We 
placed the lowest observational bound on the W fluence below 7 MeV. The relation of L- 
T, which characterizes the GRB accretion disk, is constrained. The obtained upper lim 


its are signihcantly highe r than several theoretical predictions flNagataki fc Kohril 12002 


Sekiguchi fc Shibatall201lh . However, our result is the hrst constraint that can be directly 


compared to theoretical studies. 
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DE-FG02-01ER41166, as well as other DOE grants to individual institutions, and Stichting 
Fundamenteel Onderzoek der Materie (FOM) in the Netherlands. The Kamioka Mining and 
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Fig. 1.— Schematic diagram of the KamLAND detector. The shaded region in the liqnid 
scintillator indicates the volume for the Ve analysis after the IB installation. 
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Fig. 2.— Number of long and short GRBs per three months (upper panel) and DC events 
per three months (bottom panel). The number of GRBs signihcantly increased due to the 
SWIFT satellite after December 2004. During Period I, the DC event rate is within the 
statistical uncertainty. Period II started in September 2011 and allowed for a lower energy 
threshold, increasing the number of DC events. 
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Fig. 3.— Combined background spectrum of Period I and 11 (solid red). The integral is 
7.4 X 10“^ events. The dashed curves provide our 90% upper limits under the assumption of 
a Fermi-Dirac distribution at temperature T of Fe- Each dashed curve from the left to the 
right corresponds to T = 5,10,15 MeV. 
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Fig. 4.— Fluence upper limits on Vf. from GRBs as a function of neutrino energy. Results 
from SK and SNO are also presented for comparison. SK is expected to have poorer sensi¬ 
tivity at low energies due to the detection threshold. Below 7MeV, KamLAND establishes 
the tightest limits on neutrino fluence. The slight distortion around 3-4 MeV is from the 
energy dependence of the selection efficiency (e^). 
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T [MeV] 


Fig. 5.— The KamLAND constraint on the Inminosity (L) and effective tenipera- 
tnre (T) relation in the GRB accretion disk together with several theoreti cal calcula¬ 
tions ( Nagataki fc Kohri 2002 : Sekiguchi fc Shibata 2011 : Caballero et ah 2009 ). 
























